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Abstract 

The question of naturalness is addressed in the context of gauge-mediated 
supersymmetry breaking models. Requiring that Mz arises naturally imposes 
upper limits on the right-handed selectron mass in these models that are stronger 
than in the Minimal Super symmetric Standard Model and are interesting from 
the point of view of searches at the current and future colliders. 
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How natural it is to conceive large superpartner scalar masses keeping the weak scale 
light? In this Letter we address this question of naturalness in the context of gauge-mediated 
supersymmetry breaking (GMSB) models H. There is a naturalness problem in the 
Standard Model itself, following from the assumption of the Higgs boson as a fundamental 
object, which is best cured by the elegant introduction of supersymmetry, providing a cut-off 
for the quadratic divergences associated with the Higgs. Since supersymmetry has the virtue 
of providing a radiative mechanism for spontaneous breakdown of electroweak symmmetry, 
the superpartner scalar masses cannot be arbitrarily heavy, if we require that the weak 
scale arises naturally. In other words, as those scalar masses are pushed higher, the weak 
scale still remains light only at the price of an increasingly delicate cancellation among 
the soft supersymmetry breaking masses in the theory. Thus a quantitative requirement of 
naturalness puts an upper limit on superparticle masses. In the Minimal Supersymmetric 
Standard Model with universal boundary conditions, this requirement leads to a ~ TeV scale 
upper bound on the superpartner masses §[]. We find that in GMSB models, the naturalness 
upper limits are comparatively stronger. 

In the GMSB scenario there exists a set of heavy chiral superfields, called the "messen- 
ger" fields. Supersymmetry is broken in the messenger sector due to the interaction of the 
messengers with a GUT-singlet spurion field which has a non-vanishing vacuum expectation 
value (vev) in its auxiliary component (F). The information of supersymmetry breaking is 
then transmitted to the superpartners of the quarks, leptons, gauge and Higgs bosons via the 
usual SU(3)xSU(2)xU(l) gauge interactions. The gauge-mediation attributes these models 
the virtue of automatically suppressing the Flavor-Changing Neutral Currents. In the su- 
pergravity scenario, where supersymmetry breaking in the hidden sector is communicated 
to the visible sector by gravitational strength interactions, the breaking scale y/~F is ~ 10 10 
GeV in order that the scalar masses obtained from m 2 ~ F 2 /Mp, are in the TeV range. On 
the other hand, in the GMSB models, where the scalar masses are generated radiatively, the 
analogous expressions of the squared masses have a loop-suppression factor of (a/47r) 2 and 
Mpi is replaced by the messenger mass (10 TeV&M&Mpi). Depending upon the messenger 
scale, this can lead to a smaller intrinsic supersymmetry breaking scale and consequently to 
a superlight gravitino G constituting the lightest supersymmetric particle. Of late, there has 
been a resurgence of interests [g in these models as these massless gravitinos help to provide 
a comfortable explanation to the recent CDF event with e + e~77 + missing energy in the 
final state 0. Very recently, detailed phenomenological analyses of a general class of GMSB 
models have been presented in refs. ^ and ||. 

The simplest version of the messenger models involves a set of vector-like superfields 
(Mj+Mj). These couple to a gauge singlet superfield X through the following superpotential, 



Mi and M j transform under a complete SU(5) representation and there could be n.5 copies 
of (5 + 5) and n\Q copies of (10 + 10) messengers. The quantity n = (77,5 + 3?iio) could 
at most be four from the requirement of perturbative unification. The superfield X must 
acquire non-zero vevs both in its scalar and auxiliary components through its interaction 

1 This issue has been reinvestigated with a somewhat modified criterion of naturalness in ref. Q and with 
nonuniversal boundary conditions in ref. [0. 
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with the hidden sector fields. Assuming for simplicity A = Aj (for each i), it follows from the 
superpotential written above that the messenger fermions acquire a supersymmetric massf] 
M = X(X) and the messenger scalars are split as having masses M± = M 2 ± X(Fx)- The 
effective supersymmetry breaking scale is A = (Fx) / {X) . Integrating out heavy messenger 
fields generate gaugino masses at onedoop at the M-scale, as 

Mi (M) = n ft a t (M) A g (A/M) , (2) 

where^dj = cti/Ait correspond to the gauge couplings and g(x) is the loop function. Jfc (< 1) 
arises from the fact the U(l)#-symmetry could break at a scale smaller than M rendering 
the gauginos lighter. The scalar masses arise at the twodoop level and their expressions at 
the generating scale M are given by 

3 

m 2 (M) = 2 n A 2 f (A/M) ^ C t d|(M), (3) 

i=l 

where Cj = 4/3,3/4 for the fundamental representations of SU(3), SU(2) respectively and 
zero for singlets, Q = 3Y 2 /5 for U(l) (Y = Q — T3) and f(x) is the associated loop function. 
The functions g(x) and f(x) |]l0| are very close to unity in the limit x — ► and, therefore, 
the expressions of the gaugino and scalar masses become particularly simple when A <C M, 
which is a reasonable limit that we rely upon in the rest of the papeiQ. Also, in this limit 
the gaugino and scalar masses are independent of the Yukawa couplings appearing in the 
superpotential. We do not treat n as an independent parameter, since it can always be 
absorbed in a redefinition of ^? and A as ^JnJjl — ► ^? and \fnh — > A. Actually, g and / can 
also be absorbed in a further redefinition of ^? and A. 

As has been pointed out in ref. |ll]] , there is an intrinsic "/i-problem" associated with the 
[iH u Hd term in the superpotential and the B/j>H u Hd term in the scalar potential, since these 
Peccei-Quinn symmetry-violating terms cannot be generated by gauge interactions. Both are 
generated at the same loop level in generic models rendering either fi to be at the weak scale 
and B too large beyond the naturalness bound, or B at the weak scale and /i too small to be 
phenomenologically acceptable. The problem can be cured by arranging the model in such 
a way that /j is generated at onedoop but B\x at twodoop |pd| . This keeps both fj, and B 
at the weak scale but indeed at the price of introducing additional interactions which could 
affect the Higgs boson soft mass terms. We parametrize this effect by adding two quantities 

and at the scale M in the following way: 

m 2 Hu(Hd) (M) = m 2 h {M) + □£(□2). (4) 

Soft trilinear ^4-terms are zero at M. Their radiative generation requires the simultaneous 
violation of the U(1)jj symmetry and the chiral flavor symmetry in the observable sector in 
the broken supersymmetry phase. Since the messengers do not violate observable sector 

2 We denote henceforth the superfields and their scalar components by the same symbols. 
3 Our Qi always corresponds to the one which unifies with ai and 03, i.e. it is related to the Standard 
Model ai by ai = 3«i/5. 

4 The limit A = M leads to massless messenger scalars which are definitely unwanted. Also A > M leads 
to unphysical messenger scalar masses. 
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chiral flavor, the A-terms are not generated at one-loop. They are generated only at the 
two-loop level requiring a gaugino mass insertion in the internal line ||. 

The GMSB models that we consider are characterized by the following parameters: M, 
A, ft, D 2 , d 2 ,, [x and Bq (the value of B at M). They determine the electroweak scale Mz 
and sin 2(3 through 



Tn 2 H d — f^H.. tan 2 (3 
tan 2 [3 



a fi — ILtt IU.IL U „ 

M 2 = 2—^ 5^ -~2fi 2 , (5) 



2Bu 

sin2p = — ^ -; (6) 

m% u +mj Id + 2^' 

where we have used tree-level Higgs potentials. We compute m\ using the boundary 
conditions in eqs. (||), @ and (0) and running down from M to the stop mass m^ by one- 
loop renormalization group (RG) equations^. This is enough for the purpose of studying 
naturalness. We use a^(Mz) = 0.117 and the top and bottom quark running masses at the 
top mass scale as 165 and 2.8 GeV respectively. 

Now we are all set to discuss the fine-tuning criteria. In general, if a quantity X = 
X(ai,a,2, ■■■) is determined as a function of many input parameters cij, then a measure of 
how much the parameters o&j are fine-tuned in the determination of the quantity X can be 
expressed through the parameter |3] 



A 



x 



ai dX 
X da* 



(7) 



for each cij. A measures the relative variation of X against the relative variation of ctj. In 
our case, a large fine-tuning (as defined in eq. (|7|)) corresponds to large cancellations among 
independent terms. 

Let us first examine the amount of fine-tuning that goes in the determination of sin 2(3 
through eq. ([]). For large tan/3, sin 2(3 is small and a cancellation is possible only if B = 
Bq + £? RG (where B RG is the running effect of B and that does not depend on Bq) is much 
smaller than Bq. The quantity \B RG /B\, which is in fact the same as A^| 2/3 , becomes larger 
if M increases (because in this case i? RG increases) or if tan (3 becomes larger (because in 
this case B cx sin 2(3 decreases). For example, for tan/3 ~ 50, M = 10 5 TeV, A^ 2/3 ~ 10. 

Now let us concentrate on the determination of Mz through eq. (||). Mz can be expressed 
as a function of two sets of parameters: A, M and |? in one set and fj,, d 2 and D 2 ,, which 
depend on unknown hidden sector couplings, in another. Neglecting the running of n 2 ( d ) 
(which are small), one can rewrite eq. (||) as 



fZ o rid n. u or.^ 

tan 2 /? 



M 2 =2 _tu «« --2/I 2 , (8) 



where, mi = mi, (□,, = □ ~, = 0) and 



j.2-,2 ag-^tan 2 /? 

M ~ M tan 2 /3-l ■ 1 ' 



3 This is better than running up to Mz when we consider the tree level Higgs potential. 
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Prom eq. @, it is clear that all the hidden sector couplings could effectively be expressed 
just in terms of one parameter through a redefinition of /j, 2 — > ft 2 and, therefore, there is no 
loss of generality even if we set D 2 = O 2 = in our subsequent analysis. 

Fine-tuning, in our case, corresponds to large cancellations between the two uncorrelated 
quantities in the right hand side (RHS) of eq. (||). A quantitative measure of the cancellation 
between these two terms is given by 



A = 



2/x 2 









m 



H, 



tan 2 (3 



(tan 2 (3 - 1)M§ 



1 



\M Z 
I* ' 



(10) 



Owing to our ignorance of the hidden sector couplings, A = Aff z is the best available 
measure of fine-tuning in this set-up. The rest of our paper concerns a qualitative and a 
quantitative analysis of A. 

First, we observe that for large tan/3, A ~ 2|m 2 ^/Ml|. In this case, the squared top 
Yukawa coupling (Y t ) is insensitive to tan/3, while the bottom- and tau- Yukawa couplings 
do not contribute to m/jj renormalization. Therefore, for fixed A and for tan (3 > 6, A does 
not vary with tan/3. Actually, for large tan/3, m,f R ^>me R owing to nonnegligible Y,—radiative 
corrections, so that, if we fix mf a instead of A, a small dependence of A on tan/3 creeps in. 
In any case, we hereafter confine ourselves to moderate tan (3 region. 

Apart from tan/3, A depends also on M, IfL and A (or me B ). Before demonstrating their 
exact quantitative effects, let us examine their qualitative dependences by varying them one 
at a time while keeping the others fixed. From eqs. (|2|) and @ and the RG evolution it 
follows: 



m 



u,. 



m 2 H 

flu 



m H, 



a 2 (M)+fl 2 (a 2 (M) -a 2 K))j A 2 , 

*( m t) 

dt Y t (m% + m~ + m 2 jr ) 
s v q 3 t R h u > 



(11) 
(12) 



proportional to af>, and, neglecting the Yukawa radiative corrections (just for qualitative 



where, t(Q) = 21n(M/Q) and Y t = Y t /4ir. In the RHS of eq. (H) the second term is roughly 

th 

understanding), the dependences on a 2 itself turn out to be 

'\a 2 (M) + jbf (d 2 ( 



m 



Qs,t R 



(m? 



or 



A 2 



(13) 



From eq. (|10D it follows, 



A= (J^) H{M,fi,tan[3)-l, 



(14) 



where, actually, H has a mild logarithmic dependence on A through m^, as is apparent from 
6 tan 2 /3 ^K") 



H 



tan 2 (3-1 



dt Y t 



y«i(M) + ^ 2 (a 2 -a 2 {M) 



a 2 2 {M)+fi 2 (a 2 2 {M) - a\[m L 



(15) 



Owing to the a 2 -dependence, the first term on the RHS of eq. (15) is the leading one (also, 
this term is responsible for electroweak symmetry breaking) and, consequently, H is positive. 
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The prefactor (A/Mz) 2 in eq. (|l4|) can be expressed in terms of the right-handed selectron 



mass as 

" ,2 



( a y 5 

Im z J 65? (M) 



mi R , tan 2 /3 - 1 

— — T5 sm aw ttt, r 

M| tan 2 /3 + 1 



(16) 



where, for illustration, we do not explicitly show the small U(l) RG correction. In the above 
equation, the term containing sin 2 8w corresponds to the tree-level U(l) .D-term contribution. 
For given values of tan 0, M and ft, H is determined and this could be translated into an upper 
bound of right-handed selectron mass. For example, for tan (3 = 6, M = 100 TeV and $ = 1, 
rrie R < 83 (280) GeV for A < 10 (100) (these numbers follow from exact numerical analysis). 
We emphasize that for a light selectron, the tree-level D-term relaxes the constraint. We also 
notice that the experimental lower limit on m e - fl > 45 (70) GeV from LEP1 (LEP1.5) can be 
translated through eq. (|0]) to a lower limit of A depending on M; for example, for M ~ A, 
A;> 10 (30) TeV and, on the other extreme, for M ~ M GU T, A^ 4 (15) TeV. 

In the moderate tan (3 region and for given values of M, and ft, a lower tan (3 

increases A implying a larger fine-tuning and hence a stronger constraint. Decreasing tan/3 
increases A because the leading term in H contains a factor tan 2 (3/ (tan 2 (3 — 1), Yt(mt) oc 



(tan 2 (3 + 1)/ tan 2 (3 and finally for a given mg H , (A/Mz) 2 increases (see eq. (|16|) ). 



Other parameters remaining fixed, if the U(l)# symmetry breaks at a scale lighter than 
M, |? < 1 in the minimal GMSB model (res = 1, n\Q = 0), implying a lower A and hence a 
weaker constraint. If U(l)# breaks at M, then having more messengers increases ft (from 1 
to at most 2 remaining consistent with the requirement of perturbative unification) and the 
constraint becomes stronger. 

The dependence of A on M is not apparent unlike in the cases of other parameters owing 
to several counteracting contributions. A depends in fact on M through the range of integra- 
tion, the boundary conditions on gaugino and scalar masses and the a 1 " 2 (M) dependence of 
A in eq. (]i~I| ) in conjunction with eq. (|T^). It turns out that A decreases when M — ► Mqut- 

Now we concentrate on exact numerical analysis. The dependences of A on M and ft 
are plotted in Fig. la while in Fig. lb we exhibit the limits on mg R as functions of tan/3 for 
conservative choices of M and ft extracted from Fig. la. Combining eqs. (|l|) and (H) yields 




tan 2 /3 

Mf ~ 01 " ww tan 2 p + 1 I 65? (M) 



sin 2 Q w "-^f- — } I —^——H(M, ft, tan (3) - 1 (17) 



For studying ft and M dependences, we fix tan (3 = 6. First we choose a representative value 
rrie R = 100 GeV and denote the fine-tuning parameter as Aioo- Then we plot isocurves of 
Aioo(= 3,10,20) in the log 10 (M/GeV)-ft plane. The fine-tuning for an arbitrary value of 
me R can be read from Fig. la, using 

0.21 Aioo, (18) 




which is a good approximation unless is in a physically uninteresting region, being 

much heavier than 100 GeV. It is apparent from Fig. la that in the region M ~ Mqut a 
right-handed selectron as heavy as 100 GeV does not suffer from fine-tuning problems. Also 
lowering ft reduces the fine-tuning. 
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log. (M/GeV) m- R /GeV 
(la) (lb) 

Figure 1: (a): Contour plots of Aioo (the fine-tuning parameter for nie R = lOOGeVj in 
the log 10 (M/GeV)-$ plane for tan (3 = 6; (b) naturalness upper limits on the right-handed 
selectron mass at fixed M = 10 5 TeV and ft = 1. 



In Fig. lb, we show the naturalness upper limit on mg^ as a function of tan (3 in the 
moderate tan f3 regime^ for M = 10 5 TeV (which corresponds to tHq ~ 100 eV for A = 
10 TeV g) and ft = 1. 

To conclude, we have studied the naturalness constraints in the GMSB models. The 
requirement of a natural determination of Mz in terms of the parameters of the model 
yields rather strong upper limits on the mass of the right-handed selectron. For example, for 
M = 10 5 TeV, ft = 1 and tan/3 = 6, m~ eR < 80 (230) GeV for A = 10 (100). If we choose 
A < 10 (i.e. a fine-tuning not larger than one order of magnitude), the upper limits are 
quite interesting for LEP2. We note in passing that an explanation of the CDF e + e~77 + 
missing energy event requires a cancellation of one order of magnitude in a large range of M. 
The quoted upper limits are more stringent than in the Minimal Supersymmetric Standard 
Model. In fact, in the latter scenario, the universal boundary conditions on the scalar masses 
do not permit the lightest scalar to be much lighter than the (up-type) Higgs mass. On the 
other hand, in the GMSB case, owing to the proportionality of the scalar masses to different 
gauge couplings (squared) at M, the lightest scalar becomes light enough compared to the 
Higgs mass to pose a strong fine-tuning constraint. 

We thank Riccardo Barbieri for many illuminating discussions regarding the GMSB 
models and for reading the manuscript. 

After completing our work, we have become aware of a related work in ref. [12]. 



As observed before, A is flat with respect to tan/3 for tan/3 > 6. 
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